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Abstract 

We study online scheduling problems on a single processor that can be viewed as 
extensions of the well-studied problem of minimizing total weighted flow time. In par¬ 
ticular, we provide a framework of analysis that is derived by duality properties, does 
not rely on potential functions and is applicable to a variety of scheduling problems. 
A key ingredient in our approach is bypassing the need for “black-box” rounding of 
fractional solutions, which yields improved competitive ratios. 

We begin with an interpretation of Highest-Density-First (HDF) as a primal-dual 
algorithm, and a corresponding proof that HDF is optimal for total fractional weighted 
flow time (and thus scalable for the integral objective). Building upon the salient ideas 
of the proof, we show how to apply and extend this analysis to the more general 
problem of minimizing Ylj where Wj is the job weight, Fj is the flow time 

and g is a non-decreasing cost function. Among other results, we present improved 
competitive ratios for the setting in which g is a concave function, and the setting 
of same-density jobs but general cost functions. We further apply our framework of 
analysis to online weighted completion time with general cost functions as well as 
scheduling under polyhedral constraints. 


1 Introduction 

We consider online scheduling problems in which a set of jobs JF arriving over time must 
be executed on a single processor. In particular, each job j G JF is characterized by its 
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processing time pj > 0 and its weight wj > 0, which become known after its release time 
Cj > 0. The density of j is 5j = Wj/pj, whereas, given a schedule, its completion time, Cj, 
is dehned as the hrst time t > rj such that pj units of j have been processed. The flow 
time of j is then dehned as Fj = Cj — rj, and represents the time elapsed after the release 
of job j and up to its completion. A natural optimization objective is to design schedules 
that minimize the total weighted flow time, namely the sum '^j^j processed jobs. 

A related objective is to minimize the weighted sum of completion times, as given by the 
expression '^j^jWjCj. We assume that preemption of jobs is allowed, i.e., the jobs can be 
interrupted and be resumed at a later point without any penalty. 

Im et al. [18] studied a generalization of the total weighted how time problem, in which 
jobs may incur non-linear contributions to the objective. More formally, they dehned the 
Generalized Flow Time Problem (GFP) in which the objective is to minimize the sum 
Wjg{Fj), where g : M"*" —)■ M"*" is a given non-decreasing cost function with 5 f( 0 ) = 0. 
This extension captures many interesting and natural variants of how time with real-life 
applications. Moreover, it is an appropriate formulation of the setting in which we aim 
to simultaneously optimize several objectives. We dehne the Generalized Gompletion Time 
Problem (GCP) along the same lines, with the only diherence being the objective function, 
which equals to 

A further generalization of the above problems, introduced in [5], associates each job j 
with a non-decreasing cost function gj : R"*" —)■ R’*' and gjflS) = 0; in the Job-Dependent Gen¬ 
eralized Flow Time Problem (JDGFP), the objective is to minimize the sum 'Yhj&j'^jdjiFj)- 
This problem formulation is quite powerful, and captures many natural scheduling objec¬ 
tives such as weighted tardiness. No competitive algorithms are known in the online setting; 
approximation algorithms are studied in [5l E] . 

Very recently, Im et al. [I5| introduced and studied a general scheduling problem called 
Packing Scheduling Problem (PSP). Here, at any time t, the scheduler may assign rates 
{xj{t)} to each job j G JJ. In addition, we are given a matrix B of non-negative entries. The 
goal is to minimize the total weighted flow time subject to packing constraints {Hx < 1, x > 
0}. This formulates applications in which each job j is associated with a resource-demand 
vector hj = {bij, b 2 j ,..., &mj) so that it requires an amount 6^- of the Ath resource. 

In this paper, we present a general framework based on LP-duality principles, for online 
scheduling with generalized flow time objectives. Since no online algorithm even for total 
weighted flow time is constant competitive [3] , we study the effect of resource augmentation, 
introduced by Kalyanasundaram and Pruhs (20) • More precisely, given some optimization 
objective (e.g. total flow time), an algorithm is said to be a-speed /5-competitive if it is 
/9-competitive with respect to an offline optimal scheduling algorithm of speed 1/a (here 
a > 1). 


Related work. It is well-known that the algorithm Shortest Remaining Processing Time 
(SRPT) is optimal for online total (unweighted) flow time. Becchetti et al. [B] showed that 
the natural algorithm Highest-Density-First (HDF) is (1-|-e)-speed ^^-competitive for total 
weighted flow time. At each time, HDF processes the job of highest density. 

Goncerning the online GFP, Im et al. [18] showed that HDF is (2 -|- e)-speed 0(j)- 
competitive algorithm for general non decreasing functions g. On the negative side, they 
showed that no oblivious algorithm is 0(l)-competitive with speed augmentation 2 — e, for 
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any e > 0 (an oblivious algorithm does not know the function g). In the case in which is a 
twice-differentiable, concave function, they showed that the algorithm Weighted Late Arrival 
Processor Sharing (WLAPS) is (1 + e)-speed 0(4)-competitive. For equal-density jobs and 
general cost functions [18] prove that FIFO is (1 -|- e)-speed ^-competitive. Fox et al. [10] 
studied the problem of convex cost functions in the non-clairvoyant variant providing a 
(2-|-e)-speed 0(-)-competitive algorithm; in this variant, the scheduler learns the processing 
time of a job only when the job is completed (see for example [SlElIIIlDSlEniEI]). Bansal 
and Pruhs [1] considered a special class of convex functions, namely the weighted tk norms of 
flow time, with 1 < fc < oo, and they showed that HDF is (1 -|- e)-speed 0(^)-competitive. 
Moreover, they showed how to transform this result in order to obtain an 1-speed 0(1)- 
competitive algorithm for the weighted ^k norms of completion time. 

Most of the above works rely to techniques based on amortized analysis (see also [IT] for 
a survey). More recently, techniques based on LP duality have been applied in the context 
of online scheduling for generalized flow time problems. Gupta et al. [13] gave a primal-dual 
algorithm for a class of non-linear load balancing problems. Devanur and Huang [7] used a 
duality approach for the problem of minimizing the sum of energy and weighted flow time on 
unrelated machines. Of particular relevance to our paper is the work of Antoniadis et al. [2] , 
which gives an optimal offline energy and fractional weighted flow trade-off schedule for a 
speed-scalable processor with discrete speeds, and uses an approach based on primal-dual 
properties (similar geometric interpretations arise in the context of our work, in the online 
setting). Anand et al. [1] were the hrst to propose an approach to online scheduling by 
linear/convex programming and dual htting. Nguyen [22] presented a framework based on 
Lagrangian duality for online scheduling problems beyond linear and convex programming. 
The framework is applied to several problems related to flow time such as energy plus flow 
time, £p-norm of flow time. Im et al. [15] applied dual htting in the context of PSP. For the 
weighted how time objective, they gave a non-clairvoyant algorithm that is 0(logn)-speed 
0(logn)-competitive, where n denotes the number of jobs. They also showed that for any 
constant e > 0, any (9(n^“^)-competitive algorithm requires speed augmentation compared 
to the offline optimum. 

We note that a common approach in obtaining a competitive, resource-augmented schedul¬ 
ing algorithm for how time and related problems is by hrst deriving an algorithm that is 
competitive for the fractional objective mm- An informal interpretation of the fractional 
objective is that a job contributes to the objective proportionally to the amount of its remain¬ 
ing work (see Section [2] for a formal dehnition). It is known that any a-speed /3-competitive 
algorithm for fractional GFP can be converted, in “black-box” fashion, to a (1 -|- e)a-speed 
1^/3-competitive algorithm for (integral) GFP, for 0 < e < 1 nn. Fractional objectives are 
often considered as interesting problems in their own (as in [2]). 

Contribution. We present a framework for the design and analysis of algorithms for gen¬ 
eralized how time problems that is based on primal-dual and dual-htting techniques. Our 
proofs are based on intuitive geometric interpretations of the primal/dual objectives; in par¬ 
ticular, we do not rely on potential functions. An interesting feature in our primal-dual 
approach, that dihers from previous ones, is that when a new job arrives, we may update the 
dual variables for jobs that already have been scheduled without ahecting the past portion 
(primal solution) of the schedule. Another important ingredient of our analysis consists in 
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Y.Wjg{Fj) 

Y.Wjg{Cj) 


same density 

arbitrary density 

linear 


HDF, (l + e,i^) [g 



convex 

FIFO, (l + e,l^) 

WSETF, (2-be,0(i)) HU] 
oblivious, (2 — e, P(l)) [TO] 
non-clairvoyant, (\/2 —e, 0(1)) [TU] 



concave 

differentiable 


WLAPS, (l + e,0(j^)) m 


(l + g,4(I^) 

concave 

LIFO, (1+6,1^) 

HDF, (l + e,!:^) 



general 

FIFO, (1 + e,^) m 

FIFO, (l+e,i±^) 

HDF, (2 + e,0(i)) [S] 

(7/6-e, 0(1)) m 

HDF, (1+e, 1^) 



Table 1: Summary of results for generalized flow time and completion time problems on a 
single machine. The {a, (3) notation describes algorithms which are a-speed /5-competitive. 
Our contribution (excluding our result on the PSP problem) is shown in bold. 


relating, in a direct manner, the primal integral and fractional dual objectives, without pass¬ 
ing through the fractional primal. This allows us to bypass the “black-box” transformation 
of fractional to integral solutions mm, which has been the canonical approach up to now. 
As a result, we obtain an improvement to the competitive ratio by a factor of O(-), for 
(1 -I- e)-speed. 

In Section [3] we begin with an interpretation of HDF as a primal-dual algorithm for total 
weighted flow time. Our analysis, albeit signihcantly more complicated than the known 
combinatorial one [6], yields insights about more complex problems (in the spirit of the 
dual-fitting analysis of set cover that has led to improved algorithms for metric facility 
location [I9]). Note that our approach differs from (in which the objective is to minimize 
the sum of energy and weighted flow time), even though the two settings are seemingly 
similar. More precisely, the relaxation considered in [7] consists only of covering constraints, 
whereas for minimizing weighted flow time, one has to consider both covering and packing 
constraints in the primal LP. 

In Sections 0] and [5] we expand the salient ideas behind the above analysis of HDF and 
derive a framework which is applicable to more complicated objectives. More precisely, we 
show that HDF is (1 -|- e)-speed —.competitive for GFP with concave functions, improving 
the (1 -|- e)-speed 0(^)-competitive analysis of WLAPS [H], and removing the assumption 
that g is twice-differentiable. For GFP with general cost functions and jobs of the same 
density, we show that FIFO is (1 -|- e)-speed competitive, which improves again the 
analysis in [18] by a factor of O(-) in the competitive ratio. For the special case of GFP 
with equal-density jobs and convex (resp. concave) cost functions we show that FIFO (resp. 
LIFO) are fractionally optimal, and (1 -|- e)-speed ^-competitive for the integral objective. 

In addition, we apply our framework to the following problems: i) online GGP: here, we 
show that HDF is optimal for the fractional objective, and (1 -|- e)-speed LH.Qompetitive 
for the integral one; and ii) online PSP assuming a matrix B of strictly positive elements: 
here, we derive an adaptation of HDF which we prove is 1-competitive and which reqnires 
resource augmentation max^ with Bj = maxjfejj and bj = miaiihij. 

Last, in Section [6] we extend ideas of [16], using, in addition, the Lagrangian relaxation 
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of a non-convex formulation for the online JDGFT problem. We thus obtain a non-oblivious 
(1 + e)-speed —competitive algorithm, assuming each function Qj is concave and differ¬ 

entiable. Note that this result does not rely on our framework. 

Table [1] summarizes the results of this paper in comparison to previous work. 

Notation. Let z he a job that is released at time r. For a given scheduling algorithm, we 
denote by Pr the set of pending jobs at time r (i.e., jobs released up to and including r but 
not yet completed), and by the last completion time among jobs in Pr, assuming no 
jobs are released after r. We also dehne Rr as the set of all jobs released up to and included 
r (which may or may not have been completed at r) and jTj- as the set of all jobs that have 
been completed up to time r. 

2 Linear programming relaxation 

In order to give a linear programming relaxation of GFP, we pass through the corresponding 
fractional variant. Formally, let qj{t) be the remaining processing time of job j at time t 
(in a schedule). The fractional remaining weight of j at time t is dehned as Wjqj(t)/pj. 
The fractional objective of GFP is now dehned as ~ G')- advantage of 

fractional GFP is that it admits a linear-programming formulation (in fact, the same holds 
even for the stronger problem JDGFP (see Appendix lAll). Let Xj{t) G [0,1] be a variable 
that indicates the execution rate of j G JL at time t. The primal and dual LPs are: 
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In this paper, we avoid the use of the standard transformation from fractional to integral 
GFP. However, we always consider the fractional objective as a lower bound for the integral 
one. Specihcally, we will prove the performance of an algorithm by comparing its integral 
objective to that of a feasible dual solution (D). Note that by weak duality the latter is 
upper-bounded by the optimal solution of (P), which is a lower bound of the optimum 
solution for integral GFP. 

Moreover, we will analyze algorithms that are a-speed /5-competitive. In other words, 
we compare the performance of our algorithm to an offline optimum with speed l/a (a, /9 > 
1). In turn, the cost of this offline optimum is the objective of a variant of (P) in which 
constraints ([2]) are replaced by constraints ^j(f) < 1/a for all t > 0. The corresponding 

dual is the same as (P), with the only difference that the objective is equal to ^jPj ~ 

a denote these modihed primal and dual LP’s by (Pq) and (Da), respectively. 

In order to prove that the algorithm is a-speed /5-competitive, it will then be sufficient to 
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show that there is a feasible dual solution to (-Da) for which the algorithm’s cost is at most 
/3 times the objective of the solution. 


3 A primal-dual interpretation of HDF for 

In this section we give an alternative statement of HDF as a primal-dual algorithm for 
the total weighted flow time problem. We begin with an intuitive understanding of the 
complementary slackness (CS) conditions. In particular, the primal CS condition states that 
for a given job j and time f, if Xj{t) > 0 , i.e., if the algorithm were to execute job j at time 
t, then it should be that 'j(t) = Xj — Sj{t — rj). We would like then the dual variable jit) 
to be such that we obtain some information about which job to schedule at time t. To this 
end, for any job j G JT”, we dehne the line 'jjit) = \j — 6j(t — rj), with domain [vj, oo). The 
slope of this line is equal to the negative density of the job, i.e, —6j. Our algorithm will 
always choose jit) to be equal to max{0, maxjgj-.r-<t{ 7 j(t)}} for every t > 0. We say that 
at time t the line (or the job j) is dominant if 7 j(f) = 7 (t); informally, 7 j is above y^/, for 
all j' 7 ^ j. We can thus restate the primal CS condition as a dominance condition: if a job 
j is executed at time t, then y^- must be dominant at t. 

We will consider a class of scheduling algorithms, denoted by A, that comply to the 
following rules: i) the processor is never idle if there are pending jobs; and ii) if at time r a 
new job z is released, the algorithm will hrst decide an ordering on the set Pr of all pending 
jobs at time r. Then for every f > r, it schedules all jobs in Pr according to the above 
ordering, unless a new job arrives after r. 

We now proceed to give a primal-dual algorithm in the class A (which will turn out 
to be identical to HDF). The algorithm will use the dominance condition so as to decide 
how to update the dual variables A^, and, on the primal side, which job to execute at each 
time. Note that once we dehne the A/s, the lines y/s as well as y(t) are well-dehned, as 
we emphasized earlier. In our scheme we change the primal and dual variables only upon 
arrival of a new job, say at time r. We also modify the dual variables for jobs in jZr, i.e., 
jobs that have already completed in the past (before time r) without however affecting the 
primal variables of the past, so as to comply with the online nature of the problem. 

By induction, suppose that the primal-dual algorithm A & A satishes the dominance 
condition up to time r, upon which a new job z arrives. Let qj be the remaining processing 
time of each job j G Pr at time r and \Pr\ = k. Each j G Pr has a corresponding line y^, 
once Xj is dehned. To satisfy CS conditions, each line ■jj must be dehned such that to be 
dominant for a total period of time at least Qj, in [r, cxo). The crucial observation is that, if 
a line y^ is dominant at times ti,t2, h rnust also be dominant in the entire interval [^ 1 ,^ 2 ]- 
This implies that for two jobs ji,j 2 G Pr, such that ji (resp. j 2 ) is dominant at time ti (resp. 
^ 2 ), if ti < t 2 then the slope of 7^7 must be smaller than the slope of y^j (i.e., —6j-^ < —5j^). 
We derive that A must make the same decisions as HDF. Consequently, the algorithm A 
orders the jobs in Pr in non-decreasing order of the slopes of the corresponding lines jj 
(note that the slope of the lines is independent of the A/s). For every job j G Pr, dehne 
Cj = T + ^j'l where the precedence is according to the above ordering of A. These are 

the completion times of jobs in Pr in H’s schedule, if no new jobs are released after time r; 
so we set the primal variables Xj{t) = 1 for all t G (Cj-i,Cj]. Procedure [T] formalizes the 
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Figure 1: Figure (a) depicts the situation right before r: the two lines 71 , 72 correspond to 
two pending jobs prior to the release of z. In addition, 7 (f) is the upper envelope of the two 
lines. Figure (b) illustrates the situation after the release of a third job z at time r = r^, 
the area of the shaded regions is the dual objective. In Figure (c), the area of the shaded 
regions is the primal fractional objective for the three jobs of Figure (b). 


choice of Xj for all j G Pt] intuitively, it ensures that if a job j G Pt is executed at time 
t > T then 7 j is dominant at t (see Figure [H for an illustration). 


Procedure 1 Assignment of dual variables Xj for all j G Pt 
1: Consider the jobs in Pt in increasing order of completion times Ci < C 2 < ■ ■ ■ < 
2: Choose such that 7 fc(C'fc) = 0 
3: for each pending job j = k — 1 to 1 do 
4: Choose Xj such that 'JjiCj) = 7 j+i(C'j) 


Note that the value of each Xj can be only increased after the arrival of a new job and 
the application of Procedure [TJ The following lemma shows that if no new jobs were to be 
released after time r, HDF would guarantee the dominance condition for all times t > t. 

Lemma 1 (future dominance). For Xj’s as defined by ProcedureUl and A = HDF, if job 
j G Pt is executed at time t > r, then 7 j is dominant at t. 

Proof. Suppose that the jobs in Pt are in increasing order of their completion times, i.e.. 

Cl < C 2 < ■ ■ ■ < Ck- In order to prove the lemma, we have to show that 

7j(f) < 7i+i(f) Vf > Cj (4) 

7i(f) > 7i+i(f) yte[T,Cj] (5) 

Recall that by the choice of A/s, we have that Xj — Sj(Cj — rj) = A^+i — dj^fiCj — r^+i). 

Therefore, 

~ -^i+i ~ ~ ^i+i) ^j{Cj — rj) — 6j{t — rj) 

= Xj+i — 6j^i{t — r^+i) + {6j — 6j+i){Cj — t) (6) 

Since A = HDF, we have that 5j — 5,+! >0. If f > C, then dH) follows, while if f < Cj then 

([5]) follows. □ 
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Observe that Procedure [T] modifies (increases) the Xj variables of all jobs pending at time 
r. In turn, this action may violate the dominance condition prior to r. We thus need a 
second procedure that will rectify the dominance condition for t < r. 

We consider again the jobs in Pr in increasing order of their completion times, i.e.. 
Cl < C 2 <■■■ < Ck, with k = \Pr\. We partition TZr into k disjoint sets Si, 52 ,..., Sk- 
Each set Sj is initialized with the job j G Pr, which is called the representative element of 
Sj (we use the same index to denote the set and its representative job). Informally, the set 
Sj will be constructed in such a way that it will contain all jobs a E J'r whose corresponding 
variable Aq will be increased by the same amount in the procedure. This amount is equal 
to the increase, say A^, of Xj, due to Procedure [T] for the representative job of Sj. We then 
define Procedure [2] that increases the dual variables for jobs in J'r. 


Procedure 2 Updating of dual variables Xj for all jobs j G Jr 
1: for J = 1 to A: do 
2: Add j in Sj 

3: for each job a G in decreasing order of completion times do 
4: Let b be the job such that 'ya(Ca) = 'JbiCa) 

5: Let Sj be the set that contains b 

6: Add a in Sj 

7: for each set Sj, 1 < j < /c, do 

8: Let Aj be the increase of Xj, due to Procedure [H for the representative of Sj 

9: Increase Xa by an amount of Aj for all a G Sj\ {j} 


Geometrically, the update operation is a vertical translation of the line 7 (t) for t < r. 
The following lemma shows that, if a line jj was dominant for a time t < t prior to the arrival 
of the new job at time r, then it will remain dominant after the application of Procedures [1] 
and m 

Lemma 2 (past dominance). For Xj’s as defined by both ProcedureU\ and Procedure [H and 
A = HDF, if job jEJrA Pr is executed at time t < r, then jj is dominant at t. 

Proof. The proof is based on the following three claims. 

Claim 1. Let ji and j 2 be two jobs in Pr such that Cj^ < Cj^. If A = HDF, then > Ajj. 

Proof of claim. Let z be the job released at time r, with processing time p^. Consider the 
following three cases. 

(i) Cz < Cj^ < Cj. 2 - Hence, the completion times of both ji and j 2 are delayed by pz by 
HDF. In other words, before the arrival of z the completion times of ji and j 2 were 
Cji — Pz and (5^2 — Pz, respectively. Moreover, the relative orderings of all jobs in 
Pr \ {z} (as done by the algorithm) is the same before and after the release of 2 ;. Thus, 
by Procedure [H we have that “<^ 7(^1 “Gi) ~ ~Vz — 'rjfi), where A)^ is 

the value of the dual variable of ji before the arrival of 2 ;. Hence, Aj^ = dj^pz. Similarly, 
we obtain that Aj.^ = ^j 2 Pz- Therefore, A^^ > Aj.^ as from the HDF algorithm we infer 
that dj, 






(ii) Cjj < Cj^ < Cz- In this case, Procedure [T] increases both and Xj^ by '^z{Cz — Pz — 
Tz) - lz{Cz -Tz). Therefore, = 5zPz- 

(hi) < Cz < Cj^. As in case (ii), we have A^^ = 6zPz- As in case (i), we have Aj^ = Sj^Pz- 
Therefore, Aj^ > Aj^ as from the HDF algorithm we infer that 6z > 

The claim follows. □ 

We call a set Sj critical if at least one of the following hold: Sj contains at least one job 
in Jr (i.e., Sj fl jTj- 7 ^ 0 ) or its representative job has been partially executed before time r 
(i.e., qjij) < pj). Let £ < fc be the number of critical sets. Let also Qr C Pr denote the set 
of representative jobs of critical subsets; hence \Qr\ = L Note that job z which is released 
at time r does not belong in Qr- 

Claim 2. Let ji and j 2 be two jobs in Qr such that Cj-^ < Cj^. Then > rj^. 

Proof of claim. By way of contradiction, suppose that there are two jobs ji and j 2 in Qr 
such that Cj-j < Cj.^ and < rj^. Since Cj^ < Cj^., ji has higher density than j 2 . Hence, j 2 
has not been scheduled before r because ji is active during [rj^, r] and it has higher density. 
Thus, from the dehnition of Qr, there is a job a G Sj^ H Jr for which 'ya{,Ca) = 7 ^ 2 (C'a) and 
Ca < T. Therefore, at time Ca the job j 2 is the pending job with the highest density, which 
is a contradiction as we assumed that ji is already released by time Ca and that it has higher 
density than j 2 . □ 

Let tj be the first time in which a job in the critical set Sj begins its execution for 
1 < j < Let ji,... ,ji be job indices such that tj^ < ... < tj^ < r. The following lemma 
shows a structure property of the algorithm schedule. 

Claim 3. During interval \fj^,tj^J^f), only jobs in Sj^ are executed for 1 < i < i — 1. 

Proof of claim. Consider interval \tj^,tj^J^f) and let a be the job processed at time tj^. If a 
is still pending at time r then by dehnition a is indeed job j) and Fj. consists of a singleton 
job. If a is completed at time Ca < t then by Procedure [2] there exists a job b such that 
XbiCa) = 1a{Ca) and jobs a and b belong to the same representative set. By repeating the 
same argument for job b inductively, we infer that all jobs executed in interval [tji,tj^_^_Q 
belong to the same representative set. □ 

We now continue with the proof of the lemma and we show that at each time t < r, the 
dominant job remains the same before and after the application of Procedures [1] and [2l By 
Procedure [21 the dual variables Xj of all jobs in the same critical set are all increased by 
the same amount. Moreover, by Claim [3] the jobs in the same critical set are all executed 
consecutively. Consider two jobs ji and j 2 in the same critical set Sj. Assume, without loss 
of generality, that 7 )j(t) > 7^1 (^) time t < r, where 7 '^ and 7 '^ are the lines of ji,j 2 prior 
to the arrival of z. Then 7 ji(f) = 7)1 (^) + Aj > Jj^it) + Aj = '^j^if)- Therefore, suffices to 
consider only jobs that belong in different critical sets. 

Consider the critical sets in decreasing order of completion times of their representatives, 
i.e.. Cl > C 2 > ... > Cl, and let Sj and S'^+i be any pair of consecutive critical sets. Let 
t* < T he the time at which the interval with the jobs in Sj is hnished, and a G Sj be the job 
that is executed at t*. Let also b G Sj+i be the hrst job that is executed after t*. By Claim [H 
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Aj < Aj_|_i and hence has been increased at most as mnch as A;, did (in Procedure |2]). 
Thus, 7 b is dominant right after t*. It remains to show that 7 ^ is dominant at t*. This holds 
since job b is not yet released by time t* and hence it does not affect the line 7 ^ prior to 
t*. □ 

Recall that denotes the completion time of the last pending job in (with the 

usual assumption that no job arrives after time r). The following lemma states that the 
dual variable 'y{t) has been dehned in such a way that it is zero for all t > C'max- This will 
be required in order to establish that the primal and dual solutions have the same objective 
value. 

Lemma 3 (completion). For Xj ’s defined by Procedures\I\ and\^ we have that 7 (t) = 0 for 
every t > 

Proof. From Lemmas [1] and |2l for every job j G jTj- U P^-, it holds that 7j(C'^ax) A 
By construction we have that 7(Cmax) = 0) hence 7j(Cmax) A 0. Moreover, 7 j(t) 
is non-increasing function of t, that is 7 j(t) < 7j(C'^ax) every t > Therefore, 

7 (f) = max{ 0 , maXj>.<t{ 7 j(f)}} = 0 for f > □ 

The proof of the following theorem is based on Lemmas [U [2] and [3], and it is a simplihed 
case of the proof of Theorem [2] which is given in the next section. We just note here that the 
primal and dual objectives have geometric interpretations (Figure [T](b),(c)). For each job, 
its contribution to the dual objective is the area of a trapezoid that is exactly the same as 
the contribution of the job to the primal objective. 

Theorem 1. The primal-dual algorithm A = HDF is an optimal online algorithm for the 
total fractional weighted flow time and a {1 + e)-speed ^^-competitive algorithm for the total 
(integral) weighted flow time. 

4 A framework for primal-dual algorithms 

Building on the primal-dual analysis of HDF for total weighted flow time, we can abstract 
the essential properties that we need to satisfy in order to obtain online algorithms for other 
similar problems. For the problems we consider, the primal solution is generated by an online 
primal-dual algorithm A E A which may not necessarily be HDF. In addition, each job j will 
now correspond to a curve 7 j (for the total weighted flow time problem, 7 ^ is a line), and we 
will also have a dual variable 7 (f) that will be set equal to max{ 0 ,maxjgj 7 -:rj.<i{ 7 j(f)}} for 
every f > 0. Finally, the crux is in maintaining dual variables Xj, upon release of a new job 
2 ; at time r, such that the following properties are satished: 

{VI) Future dominance. If the algorithm A executes job j at time t > t, then 7 j is 
dominant at t. 

{V2) Past dominance. If the algorithm A executes job j at time t < t, then 7 ^ remains 
dominant at t. In addition, the primal solution (i.e., the algorithm’s scheduling deci¬ 
sions) for f < r does not change due to the release of 2 ;. 
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(P3) Completion. ■y{t) = 0 for all t > 


Essentially properties (PI), (V2) and (P3) reflect that the statements of Lemmas [H [2] and [3] 
are not tied exclusively to the total weighted flow time problem. 

Theorem 2. Any algorithm that satisfies the properties {VI), {V2) and (P3) with respect to 
a feasible dual solution is an optimal online algorithm for fractional GFP and a {1 + e)-speed 
competitive algorithm for integral GFP. 

Proof. The feasibility of the solution is directly implied by the fact that Xj > 0 (since we 
only increase these dual variables) and from our dehnition of 7 (f) which implies that the 
constraints ([3]) are satished and ^{t) > 0. Let Cmax be the completion time of the last job. 
We will assume, without loss of generality, that at time t < Cmax there is at least one pending 
job in the schedule; otherwise, there are idle times in the schedule and we can apply the 
same type of analysis for jobs scheduled between consecutive idle periods. 

We will hrst show that the primal and the dual objectives are equal. Consider a job j and 
let [^ 1 ,^ 2 ], [^ 2 ,^ 3 ], • • •, [tk-i,tk] be the time intervals during which j is executed. Note that 
Xj{t) = 1 for every t in these intervals (and Xjfit) = 0 for j' 7 ^ j). Hence, the contribution of 
j to the primal (fractional) objective is 


pti+i 

'^5j / g{t-rj)dt 
i=l Vi 

By properties (PI) and (P2), the line 7 ^ is dominant during the same time intervals. Thus, 
the contribution of job j to the dual is 


^ ( rk+i ^ ( rk+i / \ ^ ( rk+i 

~ C) = / g{t - rfidt 

i=l Vi Jti \ / Jti 

since ~ Yli=i Xj{t)dt = XjPj. The hrst part of the theorem follows 

by summing over all jobs j, and by accounting for the fact that 7 (t) = 0 (from property 
(P3)), 

For the second part of the theorem, consider again the time intervals during which a job 
i is executed. The contribution of j to the integral objective is Wjg{Cj — rfi = 5jg{Cj —rfipj. 
By properties (PI) and (P2), for any t G fj+i] we have that 7 j(t) > 0. In particular, 

it holds for tk = Cj, that is Xj > Sjg{Cj —rj). Therefore, the contribution of j to the integral 
objective is 

Wjg{Cj -rj)< XjPj 

Since we consider the speed augmentation case, we will use as lower bound of the optimal 
solution the dual program that uses a smaller speed as explained in Section [2l By proper¬ 
ties (PI) and (P2), we have 7 (f) = Xj — Sjg{t — Vj) < Xj during the time intervals where the 
job i is executed. Thus, the contribution of j to the dual objective is at least 

X „ti+i ^ 

~ Z2J 'y{t)dt> XjPj - 
2=1 
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since Yl’i=i — ^jPj- From property (P3) we have that 

'y(t)dt = 0. Summing up over all jobs, the theorem follows. □ 

In what follows in this section, we apply this framework in three different problems. 

4.1 Online GCP with general cost functions 

In this section, we consider the fractional GCP and we will show that there is an optimal 
primal-dual algorithm for it. The following is a linear relaxation of the problem and its dual. 


min 



g{t)xj{t)dt 



< 1 

j€J 


Xj{t) > 0 



max 

>^jPj 

— 

poo 

/ lit)dt 





FJ 



Jo 




Vj G J 


lit) 

< 

djg{t) 

Vj 

G 

J,t> rj 



> 

0 

Vj 

G 

J 

Vt > 0 


lit) 

> 

0 

Vt 

> 

0 


Vj € J, t > 0 


We dehne 7j(t) = \j — Sg{t), based on the same arguments as in Section 3. In what 
follows, we use Procedure [3] so to define and update the values of A/s for all jobs in R^-. 


Procedure 3 Assignment and updating of Aj’s for the set Rr of all jobs released by time r. 
1: Consider jobs in R-r in increasing order of their completion times, i.e., Ci < C 2 < ■ ■ ■ < 

2: Choose Afc such that 'jkiCk) = 0 
3: for every job j = k — 1 to 1 do 
4: Let j' G Rr be the job scheduled right after Cj 

5: Choose the Xj such that 'JjiCj) = 'Jj'iCj) 


The following lemma shows that the dominance condition is always satisfied due to 
Procedure [3l 

Lemma 4. For Xj’s as defined by Procedure 0, and algorithm A = HDF, the properties 
(VI) and {V2) hold. 

Proof. Consider the jobs in Rr in increasing order of their completion times, i.e., Ci < C 2 < 
... < Ck- Let i and j -|- 1 be two consecutive jobs in this order. By the choice of dual 
variables in Procedure [3] for a time t > max{rj, rj+i} we have 

Ijif) = \ - ^j 9 if) 

= Vi - ^i+i9it) + (^i+i - ^j)i.9it) - 9 iCj)) = 7i+i(^) + (^i+i - Sfiigit) - g{Cfi) 

As we follow HDF, we have that 5j+i-5. j < 0. Since 5 ^ is a non-decreasing function, if t < Cj 
then 7 j(t) < 7^+1 (t), while if t > Cj then 7 j(t) > 7^+1 (t), and the dominance property 
follows using the same arguments as in the proof of Lemma [H 

Since the claim holds for every j, we deduce that 7 j(t) > 7 y(t) for every f G Rr and 
every time t during the execution of job j. □ 
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Property (P3) is straightforward from Procedure [3] and hence the Theorem [3] holds. 

Theorem 3. The primal-dual algorithm A = HDF is an optimal algorithm for fractional 
GCP and a (1 + e)-speed ^^-competitive algorithm for integral GCP. 

4.2 Online GFP with convex/concave cost functions and equal 
density jobs 

In this section, we consider GFP with cost function g that is either a convex or a concave 
non-decreasing function with g{0) = 0. Moreover, we assume that all jobs have the same 
density, i.e., 6j = 6 for each j G ff. For both convex and concave cases, we will show 
that there is an optimal primal-dual algorithm for minimizing the total fractional cost. For 
convex functions, this algorithm has to be the FIFO policy, whereas for concave functions 
the optimal algorithm has to be the LIFO policy. 

For both problems, we dehne jj = Xj — 6g{t — rj) (based on the same arguments as in 
Section |3] and the constraints ([3])). In what follows, we use Procedures [T] and [2] so to dehne 
and update the values of Aj’s. 

Lemma 5. For Xj’s as defined by Procedure\^ property {VI) holds if: 

(i) g is convex, all jobs have egual density and A = FIFO, 

(a) g is concave, all jobs have egual density and A = LIFO. 

Proof. We will follow the proof of Lemma [1] by showing (jl]) and (jS]). Similarly all jobs have 
the same density, we obtain 


Xj - 5g{t - rj) = Xj+i - 5g{t - r^+i) +5[{g{Cj - rj) - g{Cj - r^+i)) - {g{t - rj) -g{t - r^+i))] 

(i) Suppose that g is convex. Since Cj < Cj+i and the scheduling algorithm is FIFO, we 

have that rj < r^+i. Thus, g{Cj — rj) > g{Cj — rj+i) and g{t — rj) > g{t — Tj+i), since g 
is non-decreasing. If t > Cj, then by convexity it follows that g{Cj — rj) — g{Cj — r^+i) < 
g{t — rj) — g{t — Tj+i), and hence (jl]) holds. If t < Cj, then by convexity it holds that 

g{Cj — rj) — g{Cj — Tj+i) > g{t — rj) — g{t — rj+i), and hence (|5]) holds. 

(ii) Suppose that g is concave. Since Cj < Gj+i and the scheduling algorithm is LIFO, we 

have that rj > r^+i. Thus, g{Cj — rj) < g{Cj — rj+i) and g{t — rj) < g{t — Tj+i), since g 
is non-decreasing. If t > Cj, then by concavity it follows that g{Cj — rj) — g{Cj — r^+i) < 
g{t — rj) — g{t — r^+i), and hence (|1|) holds. If t < Cj, then by concavity it holds that 

g{Cj — rj) — g{Cj — Tj+i) > g{t — Vj) — g{t — r^+i), and hence (jS]) holds. □ 

Lemma 6. For Xj’s defined by Procedures]^ andlE, property {V2) holds if 

(i) g is convex, all jobs have egual density and A = FIFO, 

(ii) g is concave, all jobs have egual density and A = LIFO. 

Proof. We rely on the following claim. Recall that Pr denotes the set of the algorithm’s 
pending jobs at the release of z at time r. 

Claim 4. Let ji and j 2 be two jobs in P,- such that Cj^ < Cj^. Then > Ajj. 
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Proof of claim. 

(i) Suppose that g is convex and we follow the FIFO policy, then < Cj.^ < In 
this case, both 7 ^-^ and 7^2 are moved up by ^z{Cz — Pz — i^z) — lz{Cz — '^z)- Therefore, 

= ^j 2 = HgiCz - Tz) - g{Cz -Pz- Tz)). 

(ii) Suppose that g is concave and we follow the LIFO policy. Then Cz < Cj, < Cj,. 

Hence, the completion times of both ji and j 2 are delayed by p^ by the algorithm. In 
other words, before the release of 2 ; the completion times of ji and j 2 were — Pz and 
Cj 2 — Pz, respectively. Moreover, the relative ordering of all jobs in \ {z} (as done by the 
algorithm) remains the same before and after the release of 2 ;. Thus, by Procedure [H we 
have — Sg{Cj.^ — = A'^ — Sg{Cj.^ —Pz— oJ’ where A)^ is the value of the dual variable 

of ji before the release of z. Hence, = S{g{Cj.^ — — g{Cj^ — Pz — Similarly, 

^j 2 = ^( 5 '(Q 2 “’"^ 2 ) ~ 9i^j2 ~Pz ~^j 2 ))- Therefore, > Aj^ since g is concave and by the 
LIFO algorithm Oi > 02 - □ 

Note that. Claims [2] and [3] also hold for the problems we study in this section. Therefore, 
by applying the same arguments as in the proof of Lemma [2l and by replacing Claim [1] with 
Claim m we arrive at the same conclusion. □ 

The proof of the following lemma directly follows from Lemmas |5] and [ 6 ] as in the linear 
case. 

Lemma 7. For Xj’s as defined by Procedures]^ and]^ the property (VS) holds if (i) g is 
convex and all jobs have equal density; or (ii) g is concave and all jobs have equal density. 

By Lemmas [5],E] and [71 the properties (Pi), (P 2 ) and (P 3 ), respectively, are satished, and 
hence Theorem 0] holds. 

Theorem 4. The primal-dual algorithm A = FIFO (resp. A = LIFO) is an optimal online 
algorithm for fractional GFP and a {l + e)-speed ^^-competitive algorithm for integral GFP, 
when we consider convex (reps, concave) cost functions and jobs of equal density. 

4.3 PSP with positive constraint coefficients 

In this section we study the PSP problem (dehned formally in Section [1]), assuming con¬ 
straints of the form Bx < 1, and bij > 0 for every i,j- We denote by Bj the smallest and 
largest element of each column of B, respectively, i.e., bj = mini bij and Bj = maxibij. The 
primal LP relaxation of the problem, and its dual are as follows: 


poo 

jej 

)xj{t)dt 

poo 

/ Xj{t)dt>pj 

J rj 

^j^J 

Y bijXj{t) < 1 

Vt > 0, VI < i < m 

FJ 

Xj{t) > 0 

Vj € J, t > 0 


max Y, ^jPj - Y 
j&J *=i 

\j - Y^ bij^/iit) < 5j{t - rj) Mj e J,t> Vj 

i 

Xj >0 '^j £j 
7 i(t) >0 Vt > 0 


li{t)dt 
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We begin with the intuition behind the analysis, and how one can exploit the ideas of the 
analysis of HDF of Section [31 in the context of minimum total flow. The essential difference 
between the two problems is the set of packing constraints that are present in the PSP 
formulation. This difference manifests itself in the dual with the constraints Aj — bijjiit) < 
6j{t — Tj), for all j E J^,t > 0. In contrast, the dual LP of the minimum total weighted flow 
time problem has corresponding constraints Xj — 7 (f) < 6j{t — r^). At this point, one would 
be motivated to dehne /i(t) = and then view this p(f) as the variable 7 (t) in 

the dual LP formulation of the minimum total flow problem (and thus one would proceed 
as in the analysis of HDF). However, a complication arises: it is not clear how to assign the 
variables 7 j(t), for given p(f). Because of this, we follow a different way (which will guarantee 
the feasibility of the 7 j(t)’s’): Instead of satisfying the constraint Xj — J2'^i hj'^iit) < Sj(t—rj) 
we will satisfy the stronger constraint Xj — bj YlT=i Xiit) < bj{t — rj). At an intuitive level, 
we will satisfy the constraints A' — /i'(f) < 6j{t — rj), where 5' = 5j/bj and /i'(f) and A) 
are variables. Using the same scheme as in the setting of linear functions we can construct 
A' (so /i)(t)), so that the properties ("PI), (P 2 ), (P3) are satished. Particularly, if job j is 
processed at time t then p) is dominant at t, i.e., /i)(t) = A' —6j{t — rj) = = max^. 

Moreover, A) — /i'(t) < 6j(t — rj) for every job j and t > rj. 

As a last step, we need to translate the above dual variables to the PSP problem. Dehne 
Xj = X'jbj for every job j. Let j be the job processed at time t and the constraint i is a tight 
constraint at time t, i.e., bij = Bj. Set = n'it) and 7 i'(t) = 0 for i' 7 ^ i. 

Algorithm. The above discussion implies an adaptation of HDF to the PSP problem as 
follows: At any time t, we schedule job the job j which attains the highest ratio Sj/bj among 
all pending jobs, at rate Xj(t) = 1/Bj. 

Based on this algorithm, we prove Theorem [5] using similar arguments as for the proof 
of Theorem [H 

Theorem 5. For the online PSP problem with constraints Bx < 1 and bij > 0 for every i,j, 
an adaptation of HDF is ma.Xj{Bj/bj}-speed 1-competitive for fractional weighted flow time 
and maxj{(l + e)Bj/bj}-speed (1 + e)/e-competitive for integral weighted flow time. 

Proof. We hrst show the feasibility of the dual solution. For every job j and time t, we have 
Xj/bj - 6j/bj ■{t-rj)= p,'.{t) < /i'(t) = ^ 7i(^) < XI 

i i 

where the last inequality is due to bij > min* bij = bj > 0 . 

We will now bound the integral and fractional primal cost (with unit speed) by the dual 
cost (with speed bj/Bj). Consider a job j and let [^1,^2], [^2,^3], • • •, be the time 

intervals during which j is executed. Note that Xj (f) = 1/ Bj for every t during the intervals 
(and Xjflt) = 0 for f 7 ^ j). The contribution of job j to the fractional primal objective is 

m —1 pta+i m—l rta+i 

x/ bj{t-rj)bj/Bjdt 

The contribution of j to the integral primal objective is 

'^ji.Cj -rj)= Pj5j{Cj -rj)< pjXj 
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since the curve A' — — rj) > 0 for every t during which j is executed, particularly for 

t = Cj (and note that \j = bjXj and Sj = bj6j). 

Assuming that job j is processed on the machine with speed bj/Bj, the contribution of 
job j to the dual is 


h r^a + 1 m—1 pta + l / h \ 

i a=l i a=l A 3 i / 

pta+i u / \ "yA rta+i h, 

= Y = - r3)dt 

a=l A / Jta j 


where the hrst equality is due to Yl^=i = Pj! the second equality follows by Xj{t) = 

1/Bj and the dehnitions of dual variables; and the last equality holds by the dominance 
property. Hence, the contribution of job j to the fractional primal cost and that to the 
dual one (with speed bj/Bj) are equal. As the latter holds for every job, the hrst statement 
follows. 

We can bound differently the contribution of j to the dual. 


- 


bj 


rta+l 

V.Y > AjPj- 

'"j a=l dta i 


(1 + e)a 


m—1 . T 


(1 + e)B, 


E 


Xjdt 


J a=l 


= ^jPj - 


A, 


1 + e 


rta+l ^ 

j Xj{t)dt = 


where the hrst inequality due to Xj/bj = A'- > for every t during which j is 

processed; the hrst equality holds since job j is processed at rate Xj(t) = 1/Bj] and the last 
equality follows J Xj{t)dt = pj. Therefore, the contribution of job j to the integral primal 
cost is at most (1 + e)/e that to the dual objective. Again, summing over all jobs, the second 
statement holds. □ 


5 A generalized framework using dual-fitting 

In this section we relax certain properties as established in Section 0] in order to generalize our 
framework and apply it to the integral variant of more problems. Our analysis here is based 
on the dual-htting paradigm, since the analysis of Section 0] provides us with intuition about 
the geometric interpretation of the primal and dual objectives. We consider, as concrete 
applications, GFP for given cost functions g. We again associate with each job j the curve 
and set ■j(t) = ma.x{0,ma.Xj^j.,rj<t{'jj(t)}}. Then, we need to dehne how to update the 
dual variables Xj, upon release of a new job z at time r, such that the following properties 
are satished: 

(Ql) If the algorithm A schedules job j at time t > r then 'jj(t) > 0 and Xj > 'jj'it) for 
every other pending job j' at time t. 

(Q2) If the algorithm A schedules job j at time t < r, then 'jj(t) > 0 and Xj > 'jj'(t) for 
every other pending job j' at time t. In addition, the primal solution for t < r is not 
affected by the release of z. 
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(Q3) 7(^) = 0 for all t > 


Note that (Ql) is relaxed with respect to property (PI) of Section 01 since it describes 
a weaker dominance condition. Informally, (Ql) guarantees that for any time t the job 
that is scheduled at t does not have negative contributions in the dual. On the other hand, 
property (Q2) is the counterpart of (Ql), for times t < r (similar to the relation between 
(PI) and (P2)). Finally, note that even though the relaxed properties do not guarantee 
anymore the optimality for the fractional objectives, the following theorem (Theorem [6]) 
establishes exactly the same result as Theorem [2] for the integral objectives. This is because 
in the second part of the proof of Theorem 0] we only require that when j is executed at 
time t then Xj > for every other pending job j' at t, which is in fact guaranteed by 

properties (Ql) and (Q2). Therefore, the proof of the following theorem is identical with 
the one of Theorem 01 

Theorem 6. Any algorithm that satisfies the properties (Ql), (Q2) and (Q3) with respect 
to a feasible dual solution zs a (1 + e)-speed ^^-competitive algorithm for integral GFP with 
general cost functions g. 

5.1 Online GFP with general cost functions and equal-density jobs 

We will analyze the FIFO algorithm using dual htting. We will use a single procedure, 
namely Procedure 01 for the assignment of the Xj variables for each job j released by time 
r. We denote this set of jobs by Rr, and k = |P^|. 


Procedure 4 Assignment and updating of A/s for the set Rr of all jobs released by time r. 
1: Consider jobs in Rr in increasing order of completion times Ci < C 2 < ■ ■ ■ < Ck 
2: Choose Xk such that ykiCk) = 0 
3: for J = h — 1 to 1 do 

4: Choose the maximum possible Xj such that for every t>Cj, yjfi) < 7^+1 (t) 

5: if 7 j(C'j) < 0 then 

6: Choose Xj such that yj{Cj) = 0 


We will need hrst the following simple proposition: 

Proposition 1. Let i,j & J' be two jobs such that r* < Vj and suppose that there is a time 
to > p such that 7i(fo) > 7j(^o)- Then, A* > Xj. 

Proof. By the statement of the proposition we have that A* — 6g{to — r*) > Xj — 6g{to — rj). 
Since the function g is non-decreasing and r* < rj, it must be that A* > Xj. □ 

The following lemma is instrumental in establishing the desired properties. 

Lemma 8. Por every job j in Rr, Xj > 'yfiCj^i). 

Proof. Consider the jobs in Rr in increasing order of completion times. We will prove the 
lemma by considering two cases: for jobs i < j, and for jobs i > j (according to the above 
ordering). In both cases we apply induction on i. 
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Case 1: i < j. The base case i = j holds trivially. Suppose that the lemma is true 
for a job i < j, that is \j > 7 i(C'j_i). Recall that job i — 1 G i?,- is the last job that is 
completed before i. If Ai_i had not been set in line 6 of Procedure 0] then it has been set 
in line 4 of the procedure; thus, from the induction hypothesis, 7 i_i(Cj_i) < 7 i(C'i_i) < Xj. 
On the other hand, if Ai_i is set in line 6 of Procedure 0] then 7 j_i(Oj_i) = 0. Hence 
Aj > 0 = 7 j_i(C'j_i) > 7 j_i(Oj_i) where the last inequality is due to the fact that 7 i_i is a 
non-increasing function. 

Case 2: i > j. We will show the stronger claim that Xj > A, for every i > j. This suffices 
since A, > 7j(t) for every t. The base case in the induction, namely the case i = j is trivial. 

Assume that the claim is true for a job i > j, that is Xj > A*. We need to show that 
Xj > 7 i+i(C'j_i). If Xi is not set in line 6 of Procedure 0] (and thus set in line 4 instead), then 
there is time to such that 7j(to) = 7i-i-i(0) by the maximality in the choice Xj. On the other 
hand, if Aj is set in line 6 of Procedure 0] then there is at least one time ti, with Cj < ti < to 
such that 7 j(ti) > 7 j+i(ti). In both cases, from Proposition [1] we deduce that A* > Aj+i. 
Moreover, A* < Xj by the induction hypothesis. The claim follows. □ 

Lemma 9. Procedure^ satisfies properties (Ql), (Q2) and (Q3). 

Proof. From Lemma 0] it follows that properties (Ql) and (Q2) are satished. More precisely, 
observe that by line 6 of Procedure 0] we have 7 j(t) > 0 for any t < Cj. It remains to show 
that for any time t at which a job j is executed, we have Xj > ■j(t) = 'yfit). Since g 

is non-decreasing, it follows that 7 j(t) is non-increasing, hence it suffices to show the above 
only for t = Cj_i; this is indeed established by Lemma 01 

Moreover, concerning property (Q3), we observe that for all jobs j G Rr for which Xj is 
set in line 6 of Procedure 01 the property follows straightforwardly. For all remaining jobs 
in Rr, the property holds using the same arguments as in Lemma 01 for the minimum flow 
time problem. □ 

The above lemma in conjunction with Theorem [6l lead to the following result. 

Theorem 7. FIFO is a (1 -f e)-speed --competitive for integral CFP with general cost 
functions and equal-density jobs. 

5.2 Online GFP with concave cost functions 

We will analyze the HDF algorithm using dual htting. As in Section 01 we will employ two 
procedures for maintaining the dual variables Xj. The hrst one is Procedure 01 which updates 
the Aj’s for j G Ft. The second procedure updates the Aj’s for j G jTj-; this procedure is 
identical to Procedure El of Section 01 

The intuition behind Procedure Elis to ensure property (Ql) that is, 7j(t) > 0 and Xj > 
'jj'it) for all j' G Pr, which in some sense is the “hard” property to maintain. Specihcally, 
for given job j there is a set of jobs A (initialized in line 4) for which the property does not 
hold. The while loop in the procedure decreases the A values of jobs in A so as to rectify this 
situation (see line 6(ii)). However, this decrement may, in turn, invalidate this property for 
some jobs b (see line 6(i)). These jobs are then added in the set of “problematic” jobs A and 
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Procedure 5 Assignment of dual variables Xj for all j G Pr- 
1: Consider the jobs in in increasing order of completion times Ci < C 2 < ■ ■ ■ < Ck 
2: For every 1 < j <k choose Aj such that 'JjiCk) = 0 
3: for j = 2 to A: do 

4: Define A := {jobs 1 < a < j — 1 : 7 a(C'j_i) > Aj} 

5: while A / 0 do 

6: Continuously reduce Aq by the same amount for all jobs a G A until: 

(i) 3 a E A and b G Pr\A with b < a such that Aq = jbiCa-i)', then A •(— A U {6} 

(ii) 3 a € A such that 7 a(C'j_i) = Aj; then A ■(— A \ {a} 


we continue until no problematic jobs are left. One can formally argue that this procedure 
terminates. 

We will first need to show the following technical lemma: 

Lemma 10. Let a,h G LL such that 6a > 6b and suppose that there is a time to > max{ra, Vb} 
such that 7a(to) ^ 7fe(^o)- Then A^ > 7b(t) for every t > max{ra,rfc}. 

Proof. By assumption, we have Xa — 6ag{to — Va) > Xb — 6bg{to — Vb). If ra < Vb then 
g{to - Ta) > g{to - rb) so 6ag{to - ra) > 6bg{to - Vb). Hence Xa > Xb > 7b(^) for every t. 
Remains then to consider the case > r^. Since 7b(t) is non-increasing, suffices to prove 
that Xa > 7b(ra)- Since Xa> Xb — hbgifo — Vb) + hagifo — r^), it will suffice to show that 

Xb - 6bg{to - Vb) -b 6ag{to - r„) > Ab - 6bg{ra - rb) 

77 6agito - ra) > 6b{g{to - rb) - g(ra - rb)) 

Note that g is concave and 5f(0) = 0 so is sub-additive (sub-linear). Therefore, the right- 
hand side is upper bounded by 6bg{to — ra), which is at most the left-hand side. □ 

The following lemma is related to Procedure [5l which rectihes property (Ql) iteratively. 
It is not hard to see that this procedure terminates, since whenever a job a is added to 
the set A, then A^ is decreased. Note that a job may be added and removed several times 
however, A^ cannot be smaller than zero. 

Lemma 11. At the end of the for-loop for job j (line 3) of Procedure\^ 

(i) For any two jobs a,b G Pt such that 1 < a < h < j (i.e, 6a > 6b), we have Xa > Ibif) for 
all t G [Ca-i, Ca\; and 

(ii) For any two jobs a,b G Pr such that 1 <h < a < j (i.e, 6a < 6b), we have Xa > Xbif) for 
all t G [Ca-l, Ca\. 

Proof. 

(i) Proof by induction on j. For the base case {j = 1) note that in line 2 of the procedure, 
lj{Ck) = 0 for every j G Pr] hence from Lemma [TUI the base case is satished. Suppose that 
the claim holds at the end of the for-loop for job j — 1; we will call this for-loop the iteration 
for job j — 1. Consider the iteration of job j. Suppose that during this iteration, A^ has 
decreased more than Ab has, since otherwise the claim follows directly from the induction 
hypothesis. This implies that at some moment during the execution of the while loop, a G A 
and b ^ A. At that moment, we deduce that 7 b(C'j_i) < Aj < 7 a(Cj_i). However, by 
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line 6(ii), stops decreasing at the point in which Xj = 7 a(Cj_i). Therefore, at the end 
of the iteration of job j, we have < 7 a(C'j_i). Applying Lemma [TUI bv choosing 

to = Cj-i, it holds that Aq > 7fe(t) for t G [Ca_i, Cq]. 

(%%) The proof is again by induction on j. The base case, j = 1, holds since no other jobs in Pr 
have higher density than job 1. Assume that the statement holds at the end of the for-loop 
for job j — 1. During the iteration of the for-loop for job j, a set A contains jobs a < j such 
that Xj < 7 a(C'j_i). By the procedure, A^ for every a G A is decreased until Xj = 

So at the end of the iteration, Xj > 7 a(C'j_i) for every a < j. It remains to show that at 
the end of the iteration, A^ > ■jh{Ca-i) for jobs b < a < j. Let b < a < j he two arbitrary 
jobs. By the induction hypothesis, Xa > 7 f,(Ca-i) holds before the iteration. During the 
iteration, if Xa is not modihed then the inequality remains true (since A;, is not increased). 
Otherwise, a must be added to A at some moment. As Aa is decreased so probably at some 
later moment A^ = 7 b(Oa_i). However, job b will be added to A and the A-values of both jobs 
will be decreased by the same amount. Therefore, at the end of the iteration A^ > 'jb{Ca-i) 
(the induction step is done). □ 

Lemma 12. Procedures IE andlE combined satisfy properties (Q1),(Q2) and (Q3). 

Proof. We will hrst show that at the end of Procedure |5] the dual solution satishes property 
(Ql). From Lemma ITU given j E Pr, we have Xj > '^a{t) for every job a G Pt and for every 
time t G [C'j_i,C'j]. It remains to show that 7j(t) > 0 for t G [Cj-i,Cj]. Since jj is non¬ 
increasing, suffices to show that 'jjiCj) > 0. After line 2 of the procedure, 'jj{Cj) > jj^Ck) > 
0. Note that subsequently Xj may be decreased; however, we will argue that jjiCj) > 0. 
Job j may be added in A in either line 2 or in line 6(i) of the procedure, thus Xj may be 
decreased only if there is a job j' > j such that Ay < 7 j(Cj'_i). Since Cji-i > Cj and 7j(t) 
is a decreasing in t we have that Xj' < 7j(Cj). Moreover, Xj will continue decreasing until 
Ay = 7 j(Cj'_i) (line 6(ii)). Note that Xj may be decreased in many iterations and thus the 
condition Xj' = 7 j(C'y_i) may hold for more than one job j'; let j* denote the job of the last 
iteration for which Xj* = 7 j(C'y_i). Since Ay > 0 we obtain 'jjiCj) > Xj* > 0. Hence we 
showed that at the end of Procedure [5] the dual solution satishes property (Ql). 

Moreover, property (Q2) is a relaxed variant of (P2), therefore Procedure [2] guarantees 
(Q2). In addition. Procedures |5] and [2] satisfy (Q3), using very similar arguments as in the 
proof of Lemma |3l □ 

The above lemma in conjunction with Theorem [6] lead to the following result. 

Theorem 8. HDF is a (1 J- e)-speed -competitive for integral GPP with concave cost 
functions. 


6 Online JDGFP with differentiable concave cost func¬ 
tions 

We consider the online JDGFP, assuming that for each job j, the cost function pj is concave 
and differentiable. Instead of analyzing the fractional objectives and rounding to integral 
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ones (as in previous sections), we study directly the integral objective by considering a non- 
convex relaxation. Let Xj{t) be the variable indicating the execution rate of job j at time 
t. Let Cj be the variable representing the completion time of job j. We have the following 
non-convex relaxation: 


fCj 

j&j 


9 j{Cj - rj)xj{t)dt 



< 1 

J€J 


Xj(t) > 0 


yjej 

Vt > 0 

yj e j,t>o 


By associating the dual variables Xj and 7 (t) with the hrst and second constraints re¬ 
spectively, we obtain the Lagrangian dual program maxA,^ mina;^^ 7 , x, C) where the 
Lagrangian function L(A, 7 , x, C) is equal to 


L{X,-f,x,C) = 




Next, we describe the algorithm and its analysis, which is inspired by [lailE]. 

Algorithm. Let k = 1'2/e]. We write f -< j if Xji < Vj, breaking ties arbitrarily. Let 
= ^a-<j '^a.g'ait ~ ''’a) where the sum is taken over pending jobs a at time t. Note that 
as gaS are concave functions, Gj{t) is non-increasing (with respect to t). Informally, each 
term Wag'ait — r^) represents the rate of the contribution of pending job a at time t to the 
total cost. Thus Gj{t) stands for the contribution rate of pending jobs released before j 
(including j) at time t to the total cost. Whenever it is clear from context, we omit the time 
parameter. At time t, the algorithm processes job j at rate proportional to Gj{t)^ — G. 
In other words, the rate of job j at time t is = {Gj{t)^ — Gj^i{t)’‘)/G{t)’‘, where 

G{t) = J^a'^adait — Tq), and the sum is taken over pending jobs a at time t. 

Next, we define the dual variables. Dehne 7 (t) = 0. Moreover, dehne Xj such that XjPj = 

^ '^ag'aii-ra)+Wjg'j{t-rj) Yja^j dt, where Gj is the completion time 

of job j. In the following we will bound the total cost of the algorithm by the Lagrangian 
dual value. Let XF denote the total cost of the algorithm. We will use the dehnitions of the 
dual variables to derive hrst some essential properties (Lemma [13] and [TT]). 


Lemma 13. {k -I- 1) XjPj = G(t)dt = XF. 

Proof. We hrst show that {k + ^)Ylj^jPj — Consider the term Wjg'j{t — rj) in the 
sum {k -|- 1) 'Yhj ^jPj foi' f'j ^ t ^ Gj. The coefficient of this term in the sum is equal to 
^ait) + Y.a>j ^aif) = 1. Therefore, 


fCj 

{k + l)J2 \P3 = - rj)dt = Wjgj{Gj - Vj). 
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Last, note that the identity G{t)dt = ~ "^j) is straightforward from the 

dehnition of G{t). □ 

Lemma 14. For any time r > rj, it holds that Xj < Sjgj{T — rj) + Xg[t). 

Proof. The proof is by the same scheme as in [16] . We nevertheless present a complete proof, 
which follows from the following two claims, namely Claim |5] and Claim |6l 

Claim 5. For any time t > rj, 

1 / 


Pj 


Uj{t) Wag'ait - r„) + Wjg'jit - rj) Y \dt<{k + l)5jgj{T - rj) 

a<j a~ij ' 


Proof of claim. At time t, 


V. 




a-<i 


G,{tf-G,^M 

GitY 


Gj{t) < kwjg'At - rj) 


where the inequality is due to the convexity of function z^. Moreover, note that 'Yha^j — 1- 
Therefore, 

Vjit) Y '^ag'aA - To) + Wjg'j{t - rj) Y ^aif) jdt< (/c + l)wjg'j{t - rj)dt 

a^j a^j ' 


= {k + l)wjgj{T -rj) 


The claim follows. 

Claim 6. For any time r >rj, 
1 


□ 


Pj 


Y - ra) + Wjg'jit - rj) Y )dt< ( 1 + T G{t) 


a<j 


a-<j 


Proof of claim. We have 


1 

Pj 


"Ci 


Uj{t) Y - ra) + Wjg'j{t - rj) Y ^a{t) I dt 


a<j 


a^j 


< 


< 


1 

Pj Jr 
1 


(x - ra) + w.g'^t - dt 


— I w ( X - ra) + vGj+iit) ) dt 
Pj Jr 

/ 1 \ pCi 

Pj 


k 


1 + - )Cj(r) / Vj{t)dt = — ( 1 + y ]Gj{T)pj 


Pj 


k 


k 


- ( ^ + T ]Gj{.r) < ( 1 + - )G(r) 


where the hrst inequality is due to the convexity of and the second inequality follows from 
the fact that Gj{t) is non-increasing in t. Hence, the claim follows. □ 

The proof of the lemma follows by the above two claims. □ 
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Theorem 9. The algorithm is (1 + e)-speed —competitive for integral JDGFP. 
Proof. With our choice of dual variables, the Lagrangian dual objective is 


fC,- 


~ / ~ “ ^j9j{Cj - rj) dt 

- Y1 

poo 

> min J2>^3Pj- 

,• Jo , 


G{t) 


The hrst inequality is due to 7 (t) = 0 for every t and t < Cj in the integral term corresponding 
to j so gj{t — rj) < gj{Cj — rj) for every job j. The second inequality holds by Lemma HH 
In the resource augmentation model, the offline optimum has a machine of speed 
i-®-! (while the algorithm has unit-speed machine). Recall that k = [2/e]. 

Therefore, the Lagrangian dual is at least 


\Pi 

j 


1 

TT^ 





Ij^> 


k + 1 1 + ek 4(1 -|- e)2 




The algorithm has cost T so the theorem follows. 


□ 


7 Conclusion 

A promising direction for future work is to apply our framework to non-clairvoyant problems. 
It would be very interesting to obtain a primal-dual analysis of Shortest Elapsed Time First 
(SETF) which is is known to be scalable [20]; moreover, this algorithm has been analyzed 
in [10] in the context of the online GFP with convex/concave cost functions. Interestingly, 
one can use duality to argue that SETF is the non-clairvoyant counterpart of HDF; more 
precisely, one can derive SETF as a primal-dual algorithm in a similar manner as the dis¬ 
cussion of HDF in Section |3l It remains to bound the primal and dual objectives, which 
appears to be substantially harder than in the clairvoyant setting. A further open question 
is extending the results of this paper to multiple machines; here, one potentially needs to 
dehne the dual variable 7 (t) with respect to as many curves per job as machines. Last, 
we would like to further relax the conditions of the current framework in order to allow for 
algorithms that are not necessarily scalable. 

Acknowledgment. We would like to thank anonymous reviewers for pointing out an error 
in the previous version of the paper. 
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Appendix 

A LP-formulation of fractional JDGFP 

We argue that the following LP is a relaxation of JDGFP. 


min E / 9j{t - rj)xj{t)dt 
jej 


Xj(t)dt > pj 

3 

J^Xjii) < 1 

j&J 

Xj{t) > 0 


Vt > 0 


Vj e t > 0 

Consider a job j G JT”. By dehnition, the fractional weighted cost of a job j G JT” is equal 


to 


Wj{t)g'At-rj)dt = ^ / qj{t)g'Jt-rj)dt 


Pj 


= — I { 9ji^ ~ G) / Xj(u)du ) dt 

Pj Jrj 


' u=t 


By changing the order of the integrals we get 


/ \ / / \ 7 

Wj{t)gAt — Vjjdt = — 


X 


Pj 


Wi 


(u) / g'At — rj)dt du 


't=rn 


= — gj{u — rj)xj{u)du 
Pj J u=rj 


□ 
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